Parity-violating a-decay of the 3.56-MeV J^, T = 0+, 1 

state of ^Li 



Attila Csoto^ and Karlheinz Langanke^ 
^National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, 

Michigan 48824 

^ W. K. Kellogg Radiation Laboratory, California Institute of Technology, Pasadena, California 

91125 
(November 15, 1995) 



Abstract 

We study the parity-violating a + d decay of the lowest 0"'" state of ®Li in 
a microscopic three-cluster model. The initial bound and the final scatter- 
ing states are described consistently within the same model. The parity- 
violating decay width is calculated in perturbation theory using the parity- 
nonconserving nucleon-nucleon interaction of Desplanques, Donoghue, and 
Holstein (DDH). We find that the decay width is sensitive to dynamical de- 
grees of freedom which are beyond the a + p + n model, for example, a 
excitation and breakup. In our full model, which contains breathing excita- 
tions of the a particle and ^H-|-'^He rearrangement, the parity-nonconserving 
decay width is T^^^ = 3.97 • IQ-^ eV, using the DDH coupling constants. 
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I. INTRODUCTION 



Parity violation has played an important role in the understanding of the nature of the 
weak interaction. Today the weak interaction is well understood in the leptonic, semilep- 
tonic, and strange nonleptonic sectors. However, our picture of the nonstrange nonleptonic 
weak interaction, which appears in nuclear processes, e.g. in n + p scattering, is far from 
being complete The presence of the strong part of the nucleon-nucleon force makes 

the effect of the weak nucleon-nucleon interaction hardly observable. The size of the parity 
violating effect is, for example, typically in the order of 10"'' relative to the effect of the 
strong A^-A^ force. 

Despite the smallness of its effect, parity violation in a nuclear process was experimentally 
observed in 1967 by Lobashov et al. Since then, nuclear parity violation has been studied 
in p+p scattering and p + n capture, in few-nucleon systems (e.g. p + d and p + a scattering), 
in parity-mixed doublets of light nuclei, and in polarized neutron scattering on heavy nuclei. 
For two excellent review articles on the nuclear parity violation experiments, see Refs. [|T],3 . 

Complex nuclear structure and dynamics often leads to ambiguities in the theoretical 
description of nuclear parity violation involving light nuclei. However, impressive theo- 
retical progress has been achieved recently in the description of the six-nucleon systems 
PJ^J^iHiPlRCT • While realistic GFMC results for light nuclei, including for ^Li, are very 
encouraging these studies are basically restricted to bound states. However, effective 
many-body theories like the microscopic multicluster model yield a good and consistent 
description of the nuclear structure and dynamics (bound and scattering states) in the six- 
nucleon systems simultaneously This method appears thus to be well suited for the 
study of the parity-violating a + d decay of the J"^, T = 0+, 1 state of ^Li at 3.56 MeV ex- 
citation energy. As the spin-parity of the deuteron is l"*", the decay of the O"*" ^Li state into 
the a + d channel is only possible if the continuum final state is J'^ = 0~] thus this process 
violates parity. As the final state is T = 0, this process is sensitive to the isovector part of the 
parity-violating A^-A^ potential. The best experimental upper limit for the parity-violating 
a + d partial decay width is F < 6.5 ■ 10"'' eV |jll|]. The two most comprehensive theoretical 
descriptions of this process are based on the shell model |T^ and the harmonic oscillator 
cluster model [|13[ , respectively. If one assumes the same parameters for the parity- violating 
A^-A^ potential, these two studies find parity-violating decay widths which differ by orders of 
magnitude. Both of these models contain questionable approximations. In the initial O"*" 
wave function is rather schematically described by one a{pn) configuration with total spin 
and angular momentum zero. Note that the initial state has a spatially extended halo-like 
neutron-proton tail |jl4|. The ability of the shell model 0] to reproduce such a loosely 



bound state with genuine three-body {a + p + n) nature is questionable. A model which 
cannot reproduce this tail tends to compress the wave function inside, thus supposedly in- 
creasing the decay width. Moreover, both studies |T2|,|TB| use a potential model for the a + d 
scattering, which is inconsistent with the description of the bound state. 

In the present paper we study the problem in a six-body three-cluster model, which is 
practically complete in the a + p + n space. We also investigate the effects of a distortion 
and t-|-^He rearrangement on the decay width. Our model correctly describes the spatially 
extended nature of the 0"*" bound state, and uses an a + ci scattering state which is consistent 
with the bound state. We try to use a model which is as parameter-free as possible by 
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requiring the good reproduction of the properties of the subsystems (A^ + and a + N 
scattering states, channel thresholds, binding energies, etc.). 



II. MODEL 



We use the parity-nonconserving nucleon-nucleon potential of Desplanques, Donoghue, 



and Holstein |15[ , which was derived from a valence quark model. The isovector part of this 
potential reads 
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where r = |ri — r2|, fa{r) = exp{—mar) / Anr with being the vr^, p, and uj meson masses. 
Pi are the nucleon momenta, M is the nucleon mass, and cr and r are the spin and isospin 
Pauli matrices, respectively. The scalar and vector magnetic moments are xs = —0.12, and 
Xv = 3.70, and [...] and {...} denote commutators and anticommutators, respectively. We 
use the redefined coupling constants of M 



= Uq^nn/V^. Fi = -gphl/2, d = -g^hl/2. 



(2) 



The best values and reasonable ranges of the coupling constants, based on the DDH theory, 
can be found in ||l|,p^ . 



The parity- violating a + d decay width is calculated perturbatively 
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where Q{Ep) is the density of the continuum states at the final state energy Ef = (3.563 — 
1.475) = 2.088 MeV. 

Our initial ^Li state reads 
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where the indices and k denote any one of the labels a, p, and n. In (^) A is the inter- 
cluster antisymmetrizer, the $ cluster internal states are translationally invariant harmonic 
oscillator shell model states, the p vectors are the different intercluster Jacobi coordinates, 
and [...] denotes angular momentum coupling. The sum over S, li,l2, and L includes all an- 
gular momentum configurations of any significance. The last term in (^) is the t+/i=^H+^He 
rearrangement channel. The monopole breathing distortions of the a particle is considered 
by 

where the antisymmetrized ground state (z = 1) and monopole excited states (i > 1) of the 
a particle are represented by the wave functions 

$"» = ^A,,0^^., z = l,2,...,N^. (6) 

j=i 

Here (pp. is a translationally invariant shell-model wave function of the a particle with size 
parameter jSj and the parameters are to be determined by minimizing the energy of the 
a particle 0. We choose the same parameters for the wave function as in Ref. 0. We note, 
that the same model excellently reproduced the neutron halo structure of ^He 0. Thus the 
current approach is adequate to describe spatially extended systems, hke the 0"*" state of ^Li, 
which is the analog of the ^He ground state. 

The wave function of the final continuum state reads (with L = 1, = 1, and J'" = 0~), 



No. Nd 
i=2 j=2 

+ A{[[¥^%xi'{E,P,,)]^J. (7) 



Here E is the a + d relative motion energy in the center-of-mass frame. To account for 
specific distortion effects in the deuteron, we expand its wave function by 5 basis states 
{Nd = 5). The ground state ^''^ reproduces the deuteron binding energy (—2.20 MeV) and 
(point nucleon) rms-radius (1.95 fm) very well. In (^ the a + d final state is taken to be 
the time-reversed of a state with an incoming plane wave, exp(ikr), and scattered spherical 
waves. The plane wave is then projected to L = 1. The normalization of g in (|^) is chosen 
consistently with the form of the plane wave. Thus for pad oo one has 



ad) 



X (^Fi{kpad) COS 6 + Gi{kpad) sin (8) 

where k is the wave number, Fi and Gi are Coulomb functions, and 6 is the ^Pq a + d phase 
shift at energy E. 

Our model wave functions contain a large and physically most important part of the six- 
body Hilbert space. Although, our model is currently probably the closest approximation 
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to a consistent and dynamically correct full six-body description of the parity-violating 
decay process, it still has limitations. To estimate these limitations we perform a series 
of calculations in increasingly sophisticated model spaces by subsequently adding t + h 
rearrangement and a distortions to our model. 

Putting into the six-nucleon Schrodinger equation which contains a parity- 

conserving two-nucleon strong and Coulomb interaction, we arrive at an equation for the 
intercluster relative motion functions x- These functions are expanded in terms of products 
of tempered Gaussian functions exp(— 7jp^) with different ranges 7^ for each type of rela- 
tive coordinate. The expansion coefficients are determined from a variational principle. The 
scattering states are calculated from a Kohn-Hulthen variational method for the S-matrix, 
which uses square integrable basis functions matched with the correct scattering asymptotics 
16| . Then, using the resulting six-nucleon wave functions, the decay width (^) is evaluated. 



All necessary matrix elements are calculated analytically by the aid of a symbolic computer 
language P^ . 



III. RESULTS 

We consider four different model spaces with increasing level of sophistication: (i) no a 
breathing modes (A^^^ = 1) and no t + h rearrangement channel: {ai + p + n}; (ii) = 1 
and the t + h channel is included: {ai + p + n]t + h}; (iii) Na = 3 and no t + h channel: 
{as +p + n}; (iv) Na = 3 and the t + h channel is included: {a^ +p + n;t + h}. We use the 



Minnesota effective interaction |T8| as the parity-conserving strong nucleon-nucleon force. It 
has been shown ||^, that this interaction provides excellent N + N and a + N phase shifts in 
the {as + N,T + d} and {ai + A^} models by using u = 0.92 and u = 0.98 exchange mixture 
parameters of the central interaction and slightly different spin-orbit forces, respectively. 
Here we refit the u parameter in order to reproduce the experimental binding energy of 
the O"*" state (0.137 MeV relative to the a + p + n threshold) in all four models. These 
adjustments are rather small. For example, in our full model we have to change u by only 
0.7%. The weights of the {L,S) = (0,0) and (1, 1) components of the wave functions are 
between 86.5-89.5% and 13.5-10.5%, respectively in the four different models in accordance 
with the results of several different model calculations (e.g., P,p^). 



In Fig. |1| we show the 0~ phase shifts of our four models, together with the experimental 
data |]19| and the phase shift generated by the Mclntyre-Haeberli optical potential . Note 
that the two lowest-energy experimental data points are not tabulated in |T^, we read them 
off the figure of Ref. [0. All our theoretical models slightly underestimate the absolute 



value of the phase shift at the E = 2.09 MeV, which is the final energy of the decay process. 
The full model (iv) is closest to experiment. Below we investigate the sensitivity of the 
decay width on the phase shift at 2.09 MeV. 

The low-lying a + d spectrum of ^Li is reproduced well by all four models. For example, 
in the {ai+p + n,t + h} model, the J'^, T=l~^, (ground state), 3"*", 0, and 2+, ^Li states are 
at -1.421 MeV, 0.637 MeV (F = 0.012 MeV), and 4.254 MeV (F = 2.78 MeV), respectively, 
while the experimental values are -1.475 MeV, 0.71 MeV (F = 0.024 MeV), and 2.83 MeV 
(F = 1.7 MeV). All energies are relative to the a + d threshold. As in all models the deuteron 
has the correct binding energy, the position of the a + d threshold relative to the a +p + n is 



5 



always reproduced. Moreover, the ^He+p and ^Li+n thresholds are at the correct position, 
while our t + h threshold is 4-5 MeV higher than experimentally. 

In Table |I| we show the parity-violating a -\- d partial decay widths given by our four 
different model spaces, taking into account only the pion term of (|ID and using the DDH 
best value for F^. The spread among the results is almost an order of magnitude. To 
understand the origin of these big differences we performed some test calculations. We 
found that the width depends moderately on the a + d phase shift at 2.09 MeV. If we 
used a scattering wave function which reproduced the phase shift of the Mclntyre-Haeberli 
potential at 2.09 MeV, then T^^'^ was increased by roughly 50%. We checked the sensitivity 
of r on the binding energy of the 0"*" state by using the original u = 0.98 exchange mixture 
parameter in our simplest model. The binding energy decreases by 350 keV (resulting in a 
positive energy pseudo-bound state) which leads to a 15% decrease in T. Obviously, a larger 
probability of finding nucleons farther from each other decreases somewhat the magnitude 
of the internal wave function and reduces F. 

We checked the role of the orthogonal (L, S) = (0, 0) and (1,1) components in the bound 
state wave function of our four models by performing calculations for the two components 
separately. In each case the exchange mixture parameter u has been refitted to reproduce 
the correct O"*" binding energy. For the {ai + p + n} model space, the individual (0, 0) and 
(1, 1) components yield widths of 2.13 ■ 10"^" eV and 2.23 ■ 10~^ eV, respectively. This is to 
be compared with the result F = 1.06- 10^^ eV, obtained if both components are considered. 
Although we observe a strong sensitivity of F on the weight of the (L, S) = (1, 1) component, 
this cannot explain the big differences in Table | alone, as the weight of this component is 
roughly the same in all model spaces. 

In Table |I| we list the partial contributions of the (L, S) = (0, 0) and (1, 1) components 
of the bound state to the parity- violating matrix element for our four model spaces. One can 
see that the contribution of the (0, 0) component changes its sign when going from the model 
with only one a-particle basis state to the one which contains a breathing modes. This is 
obviously the consequence of some kind of cancellations taking place in the matrix element. 
These strong cancellations appear to be a warning that it is dangerous to generate bound 



and scattering states inconsistently from different models, as has been done in [12| and [13 



A similar sensitivity is found in the modelling of the beta-delayed deuteron emission process 
in ^He |21[. 

As F only moderately depends on the scattering state, the main origin of these cancella- 
tions must be in the bound state wave function. As the matrix element of the parity- violating 
interaction cannot be rewritten in a form which contains only relative motions, we cannot 
directly determine which property of the bound state wave function causes the sign change 
in the (0, 0) component. By investigating the contribution of the individual a-particle basis 
states, we find that the contributions of the excited a pseudo-states are substantial, but 



they do not cause the sign change. As in Ref. |2T|, their effect might be to shift the nodal 
positions in the bound state wave function. 

From Table II we observe that our model has a non-anticipated sensitivity to the model 
space, showing that the parity-violating decay is determined by components in the model 
space which are not well constrained by the other low-energy properties of the 6-nucleon 
systems which our models describe quite well. Nevertheless we will complete our study by 
calculating the decay width F within our most elaborate model space (iv) and using the full 
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parity-violating interaction as defined in Eqs. (l)-(2). Then 
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(9) 



In Table III we give F for the coupling constants of |T^, ||2^, and [23|, respectively. One 
observes about an order of magnitude spread among the decay widths calculated for these 
different sets of coupling constants. For comparison, Table |T| also gives F using only the 
pion contribution in the parity-violating potential. As the p and uj contributions are small, 
the decay width shows about an order of magnitude variation among the various model 
spaces. 

Finally, we compare our results with the predictions of |T^ and |T^ . In these models the 
pionic term alone gives 1.3 ■ 10~® eV and 1.1 ■ 10~^^ eV, respectively, using the DDH best 
value for the coupling constant. In |T3| only the (L, S) = (0, 0) component is included in 
the bound state wave function which might partly explain the smallness of their reported 
decay width. In |T^ the (L, S) weights are similar to ours, being 89.4% and 10.4% for (0, 0) 
and (1, 1), respectively. Compared to [0, we do not take into account the D state of the 
deuteron (our parity-conserving A^-A^ interaction is designed for a pure S-state deuteron). 
According to ||12|, inclusion of the deuteron D-state increases F by 50%. Thus our decay 
width appears to be about a factor of two smaller than the estimate given in Ref. ||12| which 
is partly due to the use of a too attractive scattering potential in [O]. In contrast to the 



present results, Ref. finds a substantial contribution of the u term in the parity- violating 
interaction to the decay width. 



IV. CONCLUSION 

In summary, we have studied the parity-violating a + d decay of the lowest 0"*" state of 
^Li within a microscopic multicluster model. We have performed a series of studies based on 
a three-cluster a + p + n model space (including all possible arrangement of the clusters and 
within each arrangement all relevant angular momenta) and its extension to additionally 
include monopole breathing modes of the a particle and the ^H-|-^He rearrangement chan- 
nel. For the parity-conserving interaction we used the Minnesota force which reproduces 
all relevant subsystem properties reasonably well. We have calculated the decay width in 
perturbation theory, using a consistent description of the 0"*" bound state and the final 0~ 
scattering state. 

We have found that the parity-nonconserving a + d decay width F^^"-^ is moderately 
sensitive to the correct reproduction of the experimental a + d phase shift and the O"*" 
binding energy. On the other hand, we have demonstrated that F^j^*^ is very sensitive to 
the presence of the 10-15% (L, S) = (1, 1) component in the bound state wave function. 
As the inclusion of a breathing modes and the ^H-|-^He rearrangement channel changes F 
considerably, we find that dynamical degrees of freedom beyond the a+p + n three-cluster 
space are quite important. The decay width in our most complete model, using the DDH 
best values for the weak coupling constants, is F^j^*^ = 3.97 ■ 10~^ eV compared to the 
experimental upper limit ^ad'" < 6.5 • 10~^ eV. 

Our most important result, however, is that the state-of-the-art microscopic multicluster 
model does not constrain the parity-violating decay width of the T = 1 state at E = 
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3.56 MeV in ^Li sufficiently well; despite the fact that this model is very successful in 
simultaneously describing other low-energy properties of the six-nucleon system. We note 
that a similar sensitivity has already been observed in cluster model studies of the beta- 
delayed deuteron emission process in ^He ||2l[]. Improved dynamical calculations in larger 
six-body model spaces would be welcome to clarify the role of other degrees of freedom 
beyond our model, for example, by including the D-state components in the deuteron and 
a-particle. Such calculations are, however, very challenging, as they have to consistently 
reproduce the spatially extended halo-like nature of the ^Li 0"*" state and the a + d scattering 
state. 
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FIGURES 

FIG. 1. ^Pq a + d phase shift as calculated in our different model spaces (dotted line: (i), 
dot-dashed: (ii), short-dashed: (iii), and solid: (iv)). The long-dashed line shows the phase shift 
as calculated from the Mclntyre-Haeberli optical potential ||2^). The dots are experimental data 
taken from Ref. |1£]. 
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TABLES 



TABLE L Parity-nonconserving a + d decay widths, calculated in the various model spaces for 
the 0"^ bound state of ^Li. Only the pion term is included in the parity-violating potential. 



Model 


Tad (eV) 


{ai +p + n} 


1.06 • 10~^ 


{ai + p + n;t + h} 


4.92 • 10^1° 


{as +p + n} 


2.15 • 10-9 


{03 + p + n;t + h} 


3.55 • 10-9 



TABLE IL Contribution of the (L, 5) = (0,0) and (1,1) components to the 
pion term of the parity- violating matrix element. The contribution is defined by 
Mi^S = -iy^27r^(E/)(^'^l^'|y/^^|^'"'^)/F^, where ^-'^i = ^-^Li + ^^Li ^^^^j ^^^^ function in 

the various model spaces. 



Model 


Moo W^) 


Mn (^) 


{ai + P + n} 


-14.5 


44.6 


{ai + p + n;t + /i} 


-12.8 


33.4 


{0:3 + p + n} 


16.0 


27.0 


{as + p + n;t + 


32.1 


23.1 



TABLE in. Parity-nonconserving a + d decay widths in our full model (iv), calculated using 
the full parity- violating potential (first column) and its pionic term only (second column). The 
calculations have been performed for various sets of weak coupling constants. 



Coupling constants 


(eV) 


^ir^^ (eV) 


DDE m 


3.97-10-9 


3.55 • 10-9 


DZ[il] 


4.28- 10-^° 


2.22 • 10-^° 


FCDH H 


1.73- 10-9 


1.21 • 10-9 
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